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ABSTRACT 

The UK Schmidt Telescope (UKST) of the Anglo- Austrahan Observatory completed 
a narrow- band Ha plus [Nil] 6548, 6584 A survey of the Southern Galactic Plane and 
Magellanic Clouds in late 2003. The survey, which was the last UKST wide- field pho- 
tographic survey, and the only one undertaken in a narrow band, is now an on-line 
digital data product of the Wide-Field Astronomy Unit of the Royal Observatory Ed- 
inburgh (ROE). The survey utilised a high specification, monolithic Ha interference 
band-pass filter of exceptional quality. In conjunction with the fine grained Tech-Pan 
film as a detector it has produced a survey with a powerful combination of area cover- 
age (4000 square degrees), resolution (~1 arcsecond) and sensitivity (<5 Rayleighs), 
reaching a depth for continuum point sources of R ~ 20.5. The main survey consists 
of 233 individual fields on a grid of centres separated by 4° at declinations below 
-1-2° and covers a swathe approximately 20° wide about the Southern Galactic Plane. 
The original survey films were scanned by the SuperCOSMOS measuring machine at 
the Royal Observatory, Edinburgh to provide the on-line digital atlas called the Su- 
perCOSMOS Ha Survey (SHS). We present the background to the survey, the key 
survey characteristics, details and examples of the data product, calibration process. 
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1 INTRODUCTION 

Ha emission from HII regions is one of the most direct opti- 
cal tracers of current star formation activity and is routinely 
used to measure star formation rates in external galaxies 
(e.g. Kennicutt 1992). In our own galaxy, HII regions are 
seen by direct UV illumination of molecular clouds from 
adjacent hot stars and as highly structured shells, bubbles 
and sheets of emission resulting from supernovae, planetary 
nebula, Wolf-Rayet stars and other stellar outflows. Some 
large-scale outflows can, in turn, be themselves a trigger of 
star formation, and their morphology is strongly influenced 
by the nature and density of the ISM into which they ex- 
pand. Ha imaging allows this to be studied in great detail in 
our immediate Galactic neighbourhood and to be detected 
at a great distance in external galaxies. The UV flux from 
hot stars also excites a more diffuse emission from the ISM, 
unconnected to current star formation and detectable over 
large areas of sky. 

The perimeters of some large emission shells appear to 
enclose the locations of more recent star formation which 
may in turn generate further supernovae and stellar winds 
(Dopita, Matthewson & Ford 1985), while their morphology 
informs the processes by which star formation is propagated 
(e.g. Elmegreen & Lada 1977, Gerola & Seiden 1978). Be- 
cause of proximity, some of these structures can present very 
large angular sizes such as Barnard's loop (probably the first 
large-scale emission structure detected in the Galaxy) which 
subtends 13° (Pickering 1890) and the Gum nebula which 
at 36°, is even larger (Gum 1952). More distant complexes 
or groups of HII regions, such as NGC 6334, can still be 
of the order 1° across and yet present fine detail on arcsec- 
ond scales (Meaburn & White 1982). Given their interaction 
with their external large-scale environment (Tenorio-Tagle 
& Palous 1987) it was clear that emission-line imaging of 
these structures required an efficient wide-field capability 
and high spatial resolution. 

On smaller scales, stellar Ha emission characterises the 
short lived, least well understood stages of stellar evolution. 



i.e. those of pre- and post-main sequence stars, planetary 
nebulae and close binary systems. Previous efforts to detect 
emission sources have either offered modest area coverage; 
e.g. the UBVI and Ha photometric surveys of Sung, Chun & 
Bessell (2000) or Keller et al. (2001) or, where a large-area 
survey has been conducted, becomes incomplete at relatively 
bright magnitudes. An example is the objective-prism sur- 
vey of Stephenson & Sanduleak (1977) which reaches only 
~ 14*'' magnitude. Such surveys are highly incomplete so 
their emission source catalogues provide only limited sam- 
ples upon which to build our understanding of these rarely 
observed phases of stellar evolution. 

Prom the above, the importance of Galactic Ha line 
emission from both stars and nebulae is evident and this has 
encouraged many surveys for HII regions in particular, e.g. 
Sharpless (1953, 1959), Gum (1955), Hase & Shajn (1955), 
Bok, Bester & Wade (1955), Johnson (1955, 1956), Rodgers, 
Campbell & Whiteoak (1960), GeorgeUn & GeorgeUn (1970) 
and Sivan (1974). These earlier surveys were limited to rel- 
atively small targeted areas or had such wide fields of view 
that small-scale detail was lost due to the low angular reso- 
lution; e.g. the survey of Sivan (1974) used 60° field diam- 
eters giving a plate scale of ~ 6° mm~^. Relatively little 
optical emission-line survey work had been done in a way 
that combined wide-angle coverage with good sensitivity and 
high resolution. These characteristics are essential to allow 
thorough examination of the morphology and interaction of 
emission regions with their environment on arcsecond to de- 
gree scales and to detect the large variety of stellar emission 
sources to suitably faint levels. 

Hence, in the mid-1990s, a number of the present au- 
thors suggested that the U.K. Schmidt Telescope (UKST) 
should be used to make a narrow band photographic Ha 
survey of the Southern Milky Way and Magellanic Clouds. 
The only previous wide-area UKST Ha material comes from 
the work by Meaburn and co-workers in the 1970s (see 
Davies, Elliot & Meaburn 1976; Meaburn 1980). They used 
a 100 A band-pass multi-element mosaic filter and fast, but 
coarse grained, 098-04 emulsion. It covered some limited ar- 
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eas close to the Galactic Plane (Meaburn & White 1982), 
but was mainly influential in the study of the ionised gas 
in the Magellanic Clouds, showing the first evidence for 
"supergiant-shells" and other largo-scale features (Davies 
et al. 1976, Meaburn 1980). There arc other recent wide- 
area Ha surveys such as that by the Virginia group in the 
northern hemisphere (VTSS - Dennison et al. 1998) and the 
Mount Stromlo group in the south (Buxton et al. 1998). Re- 
cently, Gaustad et al. (2001) have released the full "Southern 
Ha Sky Survey Atlas" (SHASSA), covering the entire south- 
ern sky. This imaging survey has rather coarse, 48-arcsecond 
pixels and strong artefacts from uncancelled stars in the 
continuum-subtracted product, but has the major benefit 
of being directly calibrated in Rayleighs. These surveys con- 
tinue the tradition of deep, low spatial resolution studies, but 
use CCDs which permit low flux densities of a few tenths of 
a Rayleigh to be achieved. 

An alternative approach, taken by the Marseille and 
Wisconsin Fabry- Perot groups in the southern and northern 
skies respectively (see Russeil et al. 1997, 1998 and Haffner 
et al. 2003 for the WHAM - Wisconsin Ha Mapper), was to 
obtain high resolution spectral (i.e. velocity) data, but again 
with low spatial resolution (e.g. 1° pixels for WHAM). 

A critical comparison between the WHAM, SHASSA 
and VTSS surveys was undertaken by Finkbeiner (2003) 
who presented a 'whole-sky' Ha map. Significantly, none 
of these major surveys offer the arcsecond spatial resolution 
of the AAO/UKST Ha survey. A summary of fundamental 
properties of these modern surveys is given in Table 1. 

2 THE AAO/UKST Ha SURVEY 

The AAO/UKST Ha survey provides a 5 Rayleigh sensi- 
tivity narrow-band survey of Galactic emission (Ha plus 
[NH] 6548, 6584A) with arcsecond spatial resolution. Hence- 
forth the survey will be refered to simply as the Ha survey 
though it is understood that this includes any [NH] emis- 
sion component that is sampled by the filter band-pass (such 
emission can completely dominate Ha for some PN types for 



example). Approximately 4000 deg^ of the Southern Milky 
Way have been covered to |6| ~ 10^ 13° together with a sep- 
arate contiguous region of 700 dcg'^ in and around the Mag- 
ellanic Clouds. Matching 3-hour Ha and 15-min broad-band 
(5900-6900A) short red (SR) exposures were taken over the 
233 distinct but overlapping fields of the Galactic Plane and 
40 fields of the Magellanic Clouds. These were done on 4° 
centres because of the circular aperture of the Ha interfer- 
ence filter which has a dielectric coating diameter of about 
305mm (~ 5.7°) deposited on a standard 356 x 356mm red 
glass (RG610) substrate (refer Section 5). The overlapping 
4° field centres enable full, contiguous coverage in Ha despite 
the circular filter aperture. Because of the slightly smaller 
effective field, a new Southern sky-grid of 1111, 4° field cen- 
tres was created (whose numbers should not be confused 
with the 893 standard 5° field centres of the UKST South- 
ern Sky Surveys). A map of the survey region in a standard 
UKST RA/DEC plot together with the new field numbers 
is available on the SHS web site^ . In the on-line UKST plate 
catalogue these fields have a 'h' prefix (e.g. hl23) to avoid 
confusion with the ESO/SERC 5° fields. 

The survey began in 1997 and took six years to com- 
plete. This latest and final UKST photographic survey was 
the first large-scale, narrow band survey undertaken on the 
telescope and is the first where the sole method of dissemi- 
nation to the community is via access to on-line digital data 
products. Preliminary survey details and results were given 
by Parker & Phillipps (1998, 2003). The present paper is 
intended as the definitive reference for the survey. We de- 
scribe the key characteristics of the survey, the on-line data 
product, some survey limitations, a flux calibration scheme, 
comparisons with other surveys and a brief overview of the 
potential for current and future scientific exploitation. 

The arcsecond resolution of the AAO/UKST Ha survey 
makes it a particularly powerful tool, not only for investi- 
gating the detailed morphology of emission features across 
the widest range of angular scales, but also as a means of 

^ http://www-wfau.roe.ac.uk/sss/halpha/ 



4 Quentin Parker et al. 

Table 1. Summajry details of various current Ha surveys 



Survey Coverage Depth Resolution Field size Filter Coverage Reference 

(sq.deg) Rayleighs (arcsec) (degrees) FWHM 

WHAMi 17000 0.15 3600 1 0.25A S > -30° Halfner et al. 2003 

SHASSA2 17000 ~ 2 48 13 x 13 32A S < 15° Gaustad et al. 2001 

VTSS^ >1000 ~2 96 10 x 10 17A 5 > -20° (if completed) Dennison et al. 1998 

SHS'' 4000 <5 1-2 5.5 x 5.5 80A 5 > +2°; |6| < 10 - 13° This paper 

IPHAS^ 1800 ~3 ~ 1 0.3 x 0.3 95A |&| < 5°; northern plane Drew et al. 2005 



Note: 1 Rayleigh = 10®/47r photons cin~^ sr^^ = 2.41 x lO"'^ erg cm~^ sr~^ at Ha 
^: http://www.astro.wisc.edu/wham; ^: http://amundsen.swaj:thmore.edu/SHASSA; 
^:http://www.phys.vt.edu/"'halpha; http://www-wfau.roe.ac.uk/ss/halpha/ 
®:http://astro.ic.ac.uk/Research/Halpha/North/ 



identifying large numbers of faint point-source Ha emitters, 
which include cataclysmic variables, T Tauri, Be and sym- 
biotic stars, compact Herbig-Haro objects and unresolved 
planetary nebulae (PNe) . Given the coincidence of the broad 
CIV/Hell blend in late-type Wolf Rayet stars, these ob- 
jects can also be detected. Most other comparative surveys 
(Table 1) are largely insensitive to point-source emitters as 
they lack spatial resolution, being optimised instead for the 
faintest levels of resolved and diffuse emission. 

On larger scales, the detailed spatial structure of the 
ionized ISM component traced by the new AAO/UKST Ha 
survey can provide key data for many studies, e.g. mapping 
of specific areas for detailed spectroscopic follow-up to ob- 
tain emission-line gas kinematics or for dynamical studies of 
star forming regions, with their implications for the energet- 
ics of the central stars. Furthermore, comparisons with other 
indicators of star formation from other wavebands should 
provide essential clues to the active mechanisms. The sur- 
vey also complements the recent Galactic Plane radio maps 
from MOST (Green et al. 1999), the new NIR maps from 
2MASS (Jarrett et al. 2000) and the mid-infrared maps from 
the MSX satellite (Price et al. 2001). 

Figure 1 presents two panels showing the 233 survey 
fields (mosaiced together by M.Read) to illustate the over- 
all survey coverage. The entire survey has been incorporated 
into an on-line mosaic within the freeware 'Zoomify' environ- 
ment (see http://www.zoomify.com) which enables prelim- 
inary survey visualisation and scanning. The lowest resolu- 
tion map can be zoomed-in to a level where each pixel repre- 



sents about 12 arcseconds. This interactive map is available 
on-line.^ This map is a factor of 18 lower in resolution than 
the full 0.67 arcsecond pixel survey data available on-line 
which should be used for serious scientific work. The success 
of this survey has led directly to a northern counterpart, 
currently underway on the 2.5m Isaac Newton Telescope on 
La Palma using a wide-field CCD camera and Ha R and I 
band imaging; the Isaac Newton telescope Photometric Ha 
Survey (IPHAS). This important survey is the subject of a 
separate paper (Drew et al. 2005) though a brief comparison 
in an overlap region is included later in Section 12. 

3 THE DETECTOR: TECHNICAL 

PANCHROMATIC FILM-BASED EMULSION 

The survey was carried out using Kodak Technical Panchro- 
matic (Tech- Pan) Estar based films (e.g. Kodak 1987). The 
superb qualities of this emulsion and its adaptation for 
UKST use has been described in detail by Parker & Mar- 
lin (1999) so only a very brief summary is given here. The 
Tech-Pan emulsion has remarkably high quantum efficiency 
for a photographic material with hypersensitised films hav- 
ing a DQE approaching 10 per cent (Phillipps & Parker 
1993). Due to its original development in connection with 
solar patrol work, it has particularly high efficiency around 
Ha. The Tech-Pan films are also extremely fine grained, 
with an inherent resolution of ~ 5/im, leading to an excel- 

^ http://surveys.roe.ac.uk/ssa/hablock/hafull.html 
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Figure 1. All 233 survey fields mosaiced together by M.Read: Top panel cover galactic longitude 1=40 to 310 degrees, bottom panel 
1=300 to 210 degrees. 



lent high-resolution imaging capability and a depth for point 
sources that exceeded that achieved for the more widely used 
glass-based IIIa-F emulsion by about a magnitude for stan- 
dard UKST R-band survey 1-hour exposures (e.g. Parker & 
Malin 1999). These factors, combined with the wide area 
coverage available to Schmidt photographic surveys, made 
Tech-Pan an ideal choice for the Southern Galactic Plane 
Ha survey. The colour term stability of Tech- Pan compared 
to the Ilia- emulsions used at the UKST is given by Mor- 
gan & Parker (2005) where these terms are shown to be 
stable, reproducible, generally small, and similar to those 
previously derived for the older Ilia- emulsions. This gives 
confidence in the survey's photometric integrity. Over the 
survey life-time, photography on a Schmidt telescope still 
offered several advantages over CCD images, especially low 
cost and very fine spatial resolution and uniformity across 
a large physical area (356 x 356 mm) giving a 40 deg^ wide 
field of view. However, a key limitation is that the detec- 
tor response is linear over only a narrow dynamic range so 
recovering and calibrating the intensity information needs 
careful treatment (see Section 11). In Figure 2 we present 



small, 3x3 arcminute regions to demonstrate the qualitative 
difference between the 3-hour Ha and 15-min SR Tech-Pan 
exposures and the standard 60-min R-band IIIarF UKST 
survey data. This region includes a newly discovered plane- 
tary nebula (PHR1706-3544) found from the Ha survey data 
as part of the Macquarie/AAO/Strasbourg Ha planetary 
nebulae catalogue (Parker et al. 2003 and 2005 in preparar 
tion). Note the improved resolution of the Tech-Pan image, 
the very similar depth of the respective exposures for point 
sources and the tighter point-spread function (psf) for the 
Tech-Pan compared to the IIIa^F emulsion. 

4 THE NARROW BAND Ha BAND-PASS 
FILTER 

To take advantage of UKST's large field of view it was nec- 
essary to obtain a physically large narrow band-pass filter to 
be placed as close as possible to the telescope's focal plane. 
The issues involved with mounting such filters with Schmidt 
telescopes has been described by Meaburn (1978) and pre- 
vious large interference filters were generally of the mosaic 
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Figure 2. 3 x 3 arcminute extracts of SuperCOSMOS data around a newly discovered PN (PHRl 706-3544) from tlie 3-liour Ha survey 
data (a - left), matching 15 minute Tech-Pan SR data (b - middle) and earlier epoch 60 minute IIIa-F R-band data (c - right). The 
new PN is only visible in the Ha image. Note the well matched depth for point sources between all three exposures and the improved 
resolution of the Tech- Pan images compared with the Hla-F exposure. 



type (e.g. Meaburn 1980). Such smaller scale interference 
filters are easier to manufacture and can be made to higher 
optical quality. However, difficulties associated with their 
mounting often load to problems of missing data in the filter 
gaps, degraded, variable resolution and lack of homogeneity 
over large survey areas, even when the optical quality of the 
elements themselves arc excellent. This was the case for the 
Meaburn mosiac filters which did not fully deliver the antic- 
ipated performance due to an unfortunate index mis-match 
in optical cement between the components which resulted in 
refiection ghosts (which can be got rid of numerically after 
scanning), coupled with the practical difficulty of mount- 
ing the components in a sandwich to eUminate optical path 
variations (Meaburn, private communication). 

Fortunately, it proved possible for the AAO to obtain 
a custom-made, exceptionally large, monolithic, thin-film 
interference filter from Barr Associates in the USA which 
avoids the problems that can be associated with mosaic 
filters. Detailed filter specifications and characteristics are 
given by Parker & Bland-Hawthorn (1998). The essential 
features are reviewed here for completeness together with 
some additional modeling of the filter profile in the converg- 
ing beam when off-axis (see Pierce 2005 for further details). 
An RG610 glass substrate was cut to 356 x 356mm (~ 6.5°), 
the standard size of UKST filters and coated with a multi- 
layer, dielectric stack to give a 3-cavity design with a clear 



aperture of ~ 305 mm diameter and with an effecive refrac- 
tive index of the equivalent monolayer of 1.34. This circular 
aperture of layered coating constitutes the interference filter 
so the corners of the square glass substrate do not behave as 
an Ha filter. Nevertheless this is probably the world's largest 
astronomical, narrow band filter. At the UKST plate scale 
this covers an on-sky area roughly 5.7° in diameter (slightly 
less than the full Schmidt field). To ensure complete and 
contiguous survey coverage with the circular aperture inter- 
ference filter it was necessary to move use 4° field centres. 

The filter central wavelength was set slightly longward 
of rest-frame Ha for two reasons, one instrumental and one 
astronomical. First, the UKST has a fast, f/2.48 converging 
beam. This leads to the interference filter 'scanning down' 
in transmitted wavelengths for oflf-axis beams compared to 
beams incident normal to the filter. Secondly we wished to 
survey positive velocity gas (in our own and nearby galax- 
ies). Given a band pass (FWHM) of 70A, we chose to centre 
the filter at 6590A in coUimated light compared to 6563A for 
rest-frame Ha. The peak filter transmission is around 90 per 
cent. Measurements of the filter at the CSIRO National Mear 
surement Laboratory in Sydney quantitatively confirmed 
the excellent conformity of the filter to our original spec- 
ifications (see Parker & Bland-Hawthorn 1998). First light 
filter images were obtained in April 1997. 
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Figure 3. The upper plot (a) shows a high resolution transmis- 
sion scan from central area of the Ha filter. The lower plot (b) 
shows a wide wavelength range transmission scan from the central 

area of the Ho filter showing the isolated narrow peak around Ha. 
Transmission at longer wavelengths is not recorded in the survey 
data as the Tech- Pan film emulsion cut-off is at 6990A and hence 
not sensitive to light at longer wavelengths. 

4.1 The filter model 

Figure 3ar-b shows two spectral scans of the filter, both taken 
near the centre using light at normal incidence. Figure 3a 
is the result of a high resolution scan around the Ha re- 
gion and shows that the bandpass is well centred on 6590A 
and has ~ 70 A FWHM. The transmission is high across the 
reasonably flat top of the bandpass, reaching over 90 per 
cent. Figure 3b is based on a scan with an extended spectral 
range from 4000 to llOOOA. The CSIRO tests show that 
the out-of-band filter transmission is 0.01 per cent or less 
up to 7600A. Figure 3b shows that the filter does transmit 
redward of 7600 A at up to ~85 per cent, but the survey 
data will be unaffected by this as the Tech-Pan film used as 
detector is insensitive beyond 6990 A. 

While this satisfies the intended performance of the fil- 



ter in light of normal incidence, in the //2.48 beam of the 
UKST, light from an object in the field centre is focused 
into a cone and enters the filter at a range of angles up to 
11.4°. The bandpass of an interference filter is blue-shifted 
for light entering at an angle. This was modelled by breaking 
down the contributions from the light cone into a series of 
concentric rings of size 1° covering the telescope beam over 
a range 0.4° to 11.4°, each entering the filter at a different 
angle. The contribution from the central part of the cone 
will not be significantly blue-shifted. The spectral shift was 
calculated for each ring according to Equation 1 adapted 
from Elliott & Meaburn (1976). 

As = Ao cos(arcsin(sin(^)//i)) (1) 

Here Ao is the chosen central wavelength for the filter 
bandpass in this case 6590 A, \g is the shifted central wave- 
length of the filter profile based on the angle, 9, of the inci- 
dent light and fi is the refractive index. A higher refractive 
index will minimise the blue shifting of the filter transmis- 
sion with incident angle of light and the filter was designed 
with this in mind. Tests performed by the CSIRO using light 
at 0°, 5° and 10° incidence found the effective refractive in- 
dex of all the layers combined, ie. the effective monolayer, is 
fj, = 1.34. This is the value used in Equation 1 to generate the 
shifts for the spectral response of the Ha filter in the UKST 
beam. These shifts are shown in Figure 4a. The solid lines 
are the shifting response curves with the most red response 
curve being applicable to light of normal incidence and the 
most blue response curve tracing the filter response to light 
entering at the most extreme angle from the telescope beam. 

In order to combine these to generate a smeared filter 
response curve which accounts for the telescope beam, each 
shifted bandpass is weighted by the area of the contributing 
ring as a fraction of the whole cone. The weighted response 
curves are shown in Figure 4b and the resulting, summed 
bandpass is shown in Figure 4c. The FWHM of this smeared 
bandpass is 80 A, centred on ~ 6550 A. 

This models the transmission of the filter in the centre 
of a survey field. Towards the edges of the field the shape 
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Figure 4. Blue-shifting response of interference filter in converg- 
ing UKST beam as one moves out from the centre of the field. 
Top plot (a) shows the shift for each concentric ring of the beam, 
the second plot (b) shows these shifted response curves weighted 
according to the area of the ring. The final plot (c) shows the 
summed, smeared out filter transmission curve. The central wave- 
length of this smeared profile is 6550 A and the FWHM is 80 A. 



of the cone changes and the maximum angle of incidence 
is over 14° which will shift the filter response further to 
the blue. The smeared out filter profile is significant as it 
permits calculation of the contribution of the contaminant 
[Nil] fines at 6548 A and 6584 A, to the flux recorded by 
the survey. Based on the smeared out filter response shown 
in Figure 4c the filter transmits HaA6563 A at 80 per cent, 



[NII]A6548 A at 82% and [NII]A6584 A at 50 per cent. Given 
that the [NIIJA6584 A line is quantum mechanically fixed to 
be three times as strong as the [NII]A6548 A line (Osterbrock 
1989), this gives a transmission of 58 per cent for any [Nil] 
emission compared with 80 per cent transmission for the Ha 
line. This is especially important when considering planetary 
nebulae (PNe) because the strength of the [Nil] lines varies 
with respect to the Ha line from PNe to PNe and will have 
a very significant impact on any calibration scheme based 
on PNe line fiux standards if not taken into account. Of 
course, for general diffuse Ha emission, the point to point 
Ha to [Nil] ratio is in general unknown without independent 
spectroscopic information, so we assume a [NIIj/Ha of 0.3, 
typically used for the warm ionised medium (e.g. Bland- 
Hawthorn et al. 1998). 



4.2 Survey depth and quality control 

The Ha films are not sky-limited after a 3-hour exposure, 
but this was chosen as a pragmatic limit which optimises 
depth, image quality and survey productivity. Field rotar 
tion and atmospheric differential refraction can adversely 
affect longer exposures (Watson 1984) which are also more 
susceptible to short-term weather and seeing variations. The 
associated 15-minute broad-band SR exposures were taken 
through the OG590 red filter. At this exposure level they are 
well matched to the depth of continuum point-sources on the 
matching Ha exposure. For completeness we include in Fig- 
ure 5 the effective SR bandpass as a function of wavelength 
obtained from a calibration spectrogram for the OG590 filter 
in combination with the Tech-Pan emulsion. 

With photographic surveys, the magnitude limit for a 
given survey field is not a fixed parameter but is a function of 
factors such as seeing, hypersensitisation and development 
of the films after exposure, emulsion batch variations and 
the brightness of the night-sky. Nevertheless, it is clear from 
comparison with the generally deeper, standard UKST R- 
band survey data, that the approximate magnitude limit for 
a typical Ha survey field in an equivalent R magnitude for 
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continuum point sources is ~ 20.5 (Arrowsmith & Parker 
2001). This value can be directly determined by examining 
the number magnitude counts from the matched Ha SR 
and R band SuperCOSMOS Image Analysis Mode (lAM) 
data (see later) for a given field and determining the point 
where completeness breaks down. As an illustration we give 
magnitude limit estimates for continuum point sources in A 
and B grade exposures of two Ha survey fields in Table 4.2. 

Additionally, the use of the same emulsion for both Ha 
and SR exposures ensures an excellent correspondence of 
their image psf's when film pairs are taken under the same 
observing conditions. The intention was to take the Ha and 
SR exposures consecutively as far possible. This greatly sim- 
plifies the inter-comparability of both types of exposure. Of 
the 233 survey fields, only 100 are in fact sequential pairs 
while most of the rest were taken a few days apart. However 
45 fields had a gap of one or more years between the Ha and 
SR survey exposures because one or other of the exposures 
had to be repeated to satisfy the stringent survey quality 
acceptance criteria. Strict quality control has been applied 
to the survey pairs by M.Hartley and S.Tritton according 
to well established criteria before any exposure is allowed to 
be incorporated into the survey. This ensures that the most 
uniform and homogeneous data set possible is created. Each 
exposure grade is determined by means of a score with '0' 
being the best and '3' being the limit for an exposure to be 
considered an 'A' grade (highest quality). The image grade 



Table 2. The depth of each of the four original images measured 
in R, equivalent magnitudes. 



Exposure 
Number 


Survey 
field 


Survey 
Grade 


Histogram num/mag peak 
(R cquiv. magnitude) 


HA18745 


h527 


A2 


20.77 ± 0.025 


HA17850 


h527 


BI3 


20.42 ± 0.025 


HA18749 


h678 


Al 


20.62 ± 0.025 


HA17935 


h678 


BIE4 


20.42 ± 0.025 



is recorded in the information and data sheets which accom- 
pany the survey data, together with a letter code to indi- 
cate which is the most significant contribution to the score. 
Long, 3-hour exposures are prone to field rotation which can 
cause image trailing (denoted by T in the image grade) , poor 
weather can lead to curtailed exposure times (U for under- 
exposed). Cosmetic defects such as emulsion faults (E), haze 
halos (H) and processing streaks (P) can also contribute to a 
poor grade. These defects can be present in either the Ha or 
the SR image. Where possible, any survey exposure which 
was not rated A grade was repeated. Unfortunately, a few 
B-grades had to be accepted into the survey though over 90 
per cent were deemed survey quality, maintaining the high 
standards set for all UKST surveys. 

5 ASTROMETRIC ACCURACY OF THE SHS 

Astrometric calibration of survey photographic material 
measured on SuperCOSMOS is discussed in Hambly et 
al. (2001). The calibration procedure consists of applying a 
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six cocfRciont (linear) plate model to measured positions of 
Tycho-2 catalogue reference stars, along with a radial distor- 
tion coefficient appropriate to Schmidt optics (i.e. tan r/r) 
and a fixed, higher order two-dimensional correction map to 
account for distortion induced by mechanical deformation 
of the photographic material when clamped in the telescope 
plate holder to fit the spherical focal surface. As demon- 
strated in Hambly et al. (2001c), this yields absolute posi- 
tional accuracy of typically ±0.2 arcseconds for glass plates. 
The SHS, on the other hand, employs film media which can- 
not be as mechanically stable as glass on the largest scales. 
However, provided a sufficiently dense grid of reference stars 
is available, it is possible to map out the unique distortion 
pattern that any one film may present. 

In order to achieve the best possible astrometry for the 
SHS, the generic SuperCOSMOS Sky Survey (SSS) astro- 
metric reduction procedure was modified by replacing the 
averaged distortion map with a correction stage where the 
individual film distortion pattern is measured with respect 
to the UCAC astrometric reference catalogue (Zacharias et 
al. 2004). In Figure 6 we show the results of comparing first- 
pass SHS astrometry (i.e. without correction of any higher 
order systematic distortion) with the UCAC catalogue for 
a single SHS film. Residuals have been averaged in 1 cm 
boxes and smoothed and filtered using a scale length of 3 
box widths. A systematic distortion pattern is clearly seen, 
and comparing with figure 1 of Hambly et al. (2001) there is 
no four-fold symmetry in the pattern, which is a characteris- 
tic of mechanical deformation of rigid glass plates. Moreover, 
similar plots for different films show different patterns, so a 
fixed correction map cannot be applied across the entire sur- 
vey film set. Figure 7(a,b) shows histograms of the residuals 
of individual UCAC standards from which Figure 6 is de- 
rived; a robustly estimated RMS (i.e. a median of absolute 
deviations scaled by 1.48, to be equivalent to a Gaussian 
sigma) is found to be about 0.4 arcseconds. Now, if the SHS 
positional data are corrected during the astrometric reduc- 
tion procedure using the map values displayed in Figure 6, 
the RMS drops to ~ 0.3 arcseconds; the new histograms of 



individual residuals are displayed in Figure 7(c,d). The value 
of ±0.3 arcseconds can be taken as indicative of the typical 
global astrometric accuracy of the SHS in either co-ordinate, 
and compares favourably with the figure quoted for the Su- 
perCOSMOS Sky Surveys (SSS) of ~ 0.2 arcseconds, given 
the higher level of crowding of the SHS fields. 

6 THE SURVEY SUPERCOSMOS DIGITAL 
DATA 

The high speed 'SuperCOSMOS' measuring machine at the 
Royal Observatory Edinburgh (e.g. Miller et al. 1992, Ham- 
bly et al. 1998) has been used to scan the Ha and SR ex- 
posure A-grade pairs at 10/ttm (0.67 arcsec) resolution. The 
same general scanning and post-processing reduction pro- 
cess is employed as for the directly analogous SuperCOS- 
MOS broad-band surveys of the Southern Sky (SSS) cur- 
rently on-line and outlined in detail by Hambly et al. (2001 
a,b,c). The user interface is broadly equivalent and the main 
features are summarised neatly in Figure 1 of Hambly et al. 
(2001a). However, due to the special nature of the survey, 
some additional processing steps and Ha specific options 
have been added to create the on-line SuperCOSMOS Ha 
Survey (SHS) described below. 

6.1 Basic characteristics of the on-line 'SHS' Ha 
Survey 

The Wide-Field Astronomy Unit (WFAU) of the Institute 
for Astronomy Edinburgh is responsible for maintaining the 
Ha survey data products. Both the Ha and SR data for 
the 233 Southern Galactic Plane survey fields are avail- 
able on-line ^. Unfortunately, there are no plans for the 40- 
field Magellanic Cloud Ha and SR survey pairs to also be 
put online. The data products are given as FITS files (see 
http://heasarc.gsfc.nasa.gov/docs/heasarc/fits.html) with 
comprehensive FITS header information detailing key pho- 
tographic, photometric, astrometric and scanning parame- 

^ http://www-wfau.roe.ac.uk/sss/halpha 
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Figure 6. Systematic astrometric distortion pattern of SHS Ha survey field h67, reduced using a standard six coeflacient linear fit plus a 
radial distortion term, when compared to the UCAC catalogue. The scale size of the vectors is 0.5 arcsec to 1 cm. Systematic positional 
errors of more than 1 arcsec (corresponding to one tick mark on either axis) are observed in the film data, e.g. in the left-hand corners. 



ters (e.g. Hambly et al. 2001b). The FITS images also have 
an accurate built-in World Co-ordinate System (WCS). This 
permits easy incorporation into other software packages such 
as the STARLINK GAIA environment for subsequent visu- 
alisation, investigation, manipulation and comparison with 
other data. The entire survey data are stored on RAID 
disks for fast access and a comprehensive set of web-based 
documentation has been provided. The pixel data map for 
each field is about 2 Gb. The scanned pixel data are pro- 
cessed through the standard SuperCOSMOS thresholded 
object detection and parameterisation software (e.g. Beard, 
MacGillivray & Thanisch 1990) to produce the associated 
Image Analysis Mode (lAM) data for each field. This pro- 
cess determines a set of 32 image-moment parameters which 
provide the astrometry, photometry and morphology of the 
detected objects. Full details of the image detection and pa- 
rameterisation are given in Hambly et al. (2001b). For the 
SHS survey, a selection of the 32 most important lAM pa- 
rameters from the merging of the Ha, SR and I band data 



for each detected image in the SHS are available and are 
given in Table 3. 

The full resolution, 10/im pixel data and associated 
lAM parameterised data for both the Ha and SR scanned 
exposures are stored on-line on a field by field basis. On the 
SuperCOSMOS web-site the scanned survey data for each 
field has the prefix 'HAL' before the survey field number 
(so Ha survey field h350 = HAL0350 for example, when re- 
ferring to the on-hne digital SuperCOSMOS data). The SR 
images have been transformed to exactly match the pixel 
grid of the master Ha exposures which permits direct im- 
age blinking and comparison between the pixel data for each 
field. The general Ha survey data products are accessed via 
a web interface that has the same look and feel as the ex- 
isting broad-band SuperCOSMOS on-line 'SSS' surveys but 
with some additional functionality. The lAM data produced 
for each field can be downloaded separately if desired or as- 
sembled into seamless catalogues on-the-fly which can cover 
several adjacent fields using the 'Get a Catalogue' option. 
The combined lAM data is organised into a full listing of 53 
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Figure 7. Histograms of residuals between SHS and UCAC astrometry (see text) for: (a,b) uncorrected positions, and (c,d) positions 
corrected during the astrometric reduction procedure using the distortion map shown in Figure 6. The accuracy is quantified by a 
robustly estimated RMS residual, in cither co-ordinate, and shows a ~ 30% improvement when the correction map is employed. Prom 
this analysis, SHS absolute astrometry is typically accurate to ~ 0.3 arcsec. 



image parameters or a more manageable subset of the most 
useful 32 as in Table 3. A set of 'expert' options are also 
available to further select catalogue extraction parameters. 
A special feature to create a difference image of each field fol- 
lowing variable image psf matching techniques developed by 
Bond et al. (2001) also exists to permit large-scale resolved 
emission maps to be created with reduced artefacts from un- 
cancelled stars. This can be computationally intensive and 
so is not generally available without prior arrangement with 
the Wide-Field Astronomy Unit. For most applications sim- 
ple quotient imaging between the Ha and SR pixel data is 
sufiicient due to the well matched psf 's and depth. 

A 16x blocked-down version of each field is also avail- 



able as both a GIF image and as a FITS file which has the 
WCS built in to the FITS header. These whole field maps 
can be studied to select smaller regions of interest for ex- 
traction at full resolution using the 'Get an Image' option. 
The full resolution pixel data access limit is currently set at 
9000 arcmin^ with regions downloaded as FITS files (also 
with WCS) and both the SR and Ha data for the same re- 
gion can be downloaded simultaneously. Areas for extraction 
can be chosen via equatorial (J2000 or B1950) or galactic 
(l,b) co-ordinates in a variable mx n arcminute rectangular 
region format. A clickable map of the current fields on-line 
enables individual field details to be displayed prior to view- 
ing the blocked full field image. A batch mode enables large 
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Table 3. The selected 32 lAM parameters used in the merged SHS catalogue data 



Number 


Name 


Type 


Description 


Units 


\ 


RA 






1 duidiio I J. J. _L o J 


2 


DEC 




Celestial Declination 


radians (^i'^lTS^ 


3 


EPOCH 


Real 


Epoch 


year 


4 


12 


Real 


lip 1 nr'i'ir' 1(~iti cti Hf> 

vjr<xiaj\_> LiH_> iwiitiiLrU.L_is_. 


Llt^^ , Li til CCD 


5 


b2 


Real 


Galactic latitude 


Hop Hpptpp^ 


R 

o 


r\,_na 


Real 


Hoi ecjuivalcnt R-iiiag 


mags 




SR2 


Real 


YY\ Oli'^'V\^Y^ (T ^ l-< yy'\ c\ (t\^ ^ i' ^'\ t'\ I ^ l-s v i 
lllctH_-illlig OlV. illcL^lil L U.U(_- lOlV.^) 


mags 


g 


SRI 


Real 






g 


I 


Real 


nivct pnr^r*n T TnncrnTl'nriP 

J. 11 oil K!\J\J\jlL 1. lllCLgllllj uuc? 


mags 


10 


AREAJia 


Tntpp'PT 


TVit a 1 a rpa TT/T" i m a o'p 

-L yj uciii cti. ^cu J. J.t_^ iiiidr^^ 


ni vpI 


11 


AREA_SR2 


Tntpcpr 


TVita] arpfi SR*? iTnfiP"p 


ni ypI *5 


12 


AREA_SR1 


Integer 


Total area SRI image 


pixels 


13 


ARRA T 


±11 LL>ti^i 


_LW L cxi <ni \_.<n 1 Illl<Xti^ 


nivplc; 

^lAClO 


14 


ET.T, Ha 


Real 


1 1 (Xl IJllOj IlllcXti^ L,llll_/LH_>ILry 




15 


ELL_SR2 


Real 


SR2 image elliptic! ty 




16 


ELL_SR1 


Real 


SRI image ellipticity 




17 


ELLJ 


Real 


T imap'p pllintifitv 

A iiiiCL^^ ^iiiu iji^i yj\ 




18 


PRFSTT_Ha 


Real 


HcK N(0,l) profile classification statistic 


0.001 sigma 


19 


PRFSTT SR2 


Real 


STR ^ IV f nil 'nrrtfilp plfm*?ififatinn *?tati^tip 

kJ-Llj^ l^\V^^-Li 1-^ Wlll^ ^ICtoOllll^CliUl^ll O vClrvlO un^ 


Oni siema 

\J m\J\J ±. Olg^lllClr 


20 


PRFSTT.SR1 


Real 


SRI N(0,1) profile classification statistic 


0.001 Sigma 


21 


PRFSTTJ 


Real 


I N(0,1) profile classification statistic 


0.001 sigma 


22 


BLEND.Ha 


Integer 


Ha Dcblcnding flag (0 if not dcblcnded) 




23 


BLEND.SR2 


Integer 


SR2 Dcblcnding flag (0 if not dcblcnded) 




24 


BLEND.SR1 


Integer 


SRI Dcblcnding flag (0 if not deblended) 




25 


BLENDJ 


Integer 


I Dcblcnding flag (0 if not deblended) 




26 


QUAL.Ha 


Integer 


Ha image Quality flag 




27 


QUAL.SR2 


Integer 


SR2 image Quality flag 




28 


QUAL.SR1 


Integer 


SRI image Quality flag 




29 


QUALJ 


Integer 


I image Quality flag 




30 


PA 


Integer 


Celestial position angle (selected band) 


degrees 


31 


CLASS 


Integer 


Classiflcation flag (selected band) 




32 


FIELD 


Integer 


SHS field number (4° centres) 





numbers of thumb-nail images to be extracted around ob- 
jects of interest with the option to return Ha and/or SR 

postscript plots of the extracted images. An option to ap- 
ply a 'Flat-Field' to the Ha pixel data in intensity space is 
included to permit correction of the non-uniformities in the 
me£isured exposures arising from the excellent but slightly 
varying Ha filter transmission profile. This has been shown 
to work effectively and is described in Section 9. A radius 
of 153 mm (~ 2.85 degrees) from each survey field centre 
has boon adopted as the region with good data (< 15 per 
cent correction factors). The 'good data' radius from each 
scanned Ha field centre has been used in creating a confi- 
dence map which is incorporated into the extracted FITS 
image as an additional FITS extension (extension [3], e.g 
test-image. fits[3]). This can be used to flag areas of the ex- 
tracted image that might not be quite as good as others. 



Currently this has values of 100 for regions extracted inte- 
rior to this radius and for regions outside. 

6.2 Incorporation of the SSS '/' band data 

The I AM catalogue downloaded directly via the 'Get a Cat- 
alogue' option or as incorporated in the FITS table exten- 
sion to the dowloaded pixel data via the 'Get an Image' 
option, contains information not just from the Ha and the 
SR images but also from the SERC-I (near infrared) survey 
which has been carefully matched in with the SHS data. The 
I-band data is particularly useful when searching for point- 
source Ha emitters as it can help to eliminate contamination 
from late-type stars. However there are issues that any user 
should be awaxe of when combining the I-band data with the 
SHS magnitudes. When calibrating UKST data, positional 
and magnitude-dependent systematic errors axe present as a 
result of variations in emulsion sensitivity and vignetting to- 
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wards the image corners (Harnbly et al. 2001b). The I-band 
survey was taken on standard ESO/SERC 5° field centres, 
so the vignetting will have a diflterent effect on the two sets 
of photometry at a given survey location. Furthermore, the 
I-band data are calibrated to relatively few standards. These 
differences axe evident when looking at a plot of SR mag- 
nitude versus R — I colour derived from the SR and paired 
I-band photometry. Figure 8 shows two colour magnitude 
diagrams (CMD) for stars taken from a 10 arcminute re- 
gion centred on the middle of SHS field hi 109 in Mono- 
ceros where the low galactic reddening of E{B — V) =0.24 
(Schlegel, Finkbeiner & Davis 1998) should leave the R — I 
colour roughly constant for much of the observed magni- 
tude range. Figure 8a shows the raw result, where a large, 
unphysical variation of 3.5 magnitudes is seen in the R — I 
colour from the brightest to the faintest objects. This is re- 
moved as a first-order correction from the survey data by 
selecting a master colour, in this case the SR, and correct- 
ing the I-band across all survey fields. Note however that 
at fainter limits one in fact expects redder R — I colours 
as such stars are likely to be further away and prone to 
be more dust reddened or intrinsically fainter and therefore 
more likely to be late types. Hence some modest slope is 
expected. Figure 8b shows the CMD for the same patch of 
sky after the colour correction has been applied, giving a re- 
sult in better agreement with expectation. The data can be 
downloaded in corrected or uncorrected form via an option 
in the "Expert" parameters of the SHS website. It is im- 
portant to ensure that the I-band correction is not applied 
inappropriately, i.e. in a field of intrinsically high reddening, 
as such a correction will remove genuine features from the 
data. It should be used with caution. 

7 SHS POINT SOURCE PHOTOMETRY 

A significant advantage of the SHS data over its rivals is the 
ability to detect point sources which have been photomet- 
rically calibrated to CCD standards (e.g. Boyle et al. 1995, 
Croom et al 1999). With measurements of isophotal mag- 
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Figure 8. Plots of SR magnitude versus R — I, from a 10 ar- 
cminute square region in Monoceros, before colour correction for 
positional and magnitude dependent errors (a - upper ) and after 
(b - lower). 

nitude and object classifications, it is possible to apply a 
photometric calibration to the Ha and SR films by compar- 
ing the SuperCOSMOS raw magnitudes of stars from the 
Tycho-2 Catalogue (Hog et al. 2000) and the Guide Star 
Photometric Catalogue (Lasker et al. 1988). These in turn 
are checked against photometric standards derived from the 
CCD observations given by Croom et al. (1999) and Boyle 
et al. (1995). The narrow-band Ha images are calibrated to 
an 'R-equivalent' scale. The 3 hour Ha and 15 minute SR 
exposures are matched so that both reach similar depths of 
thr ^ 20.5 for point sources. Where an object is detected in 
one band but not in the other a default value of 99.999 is 
given in the catalogue data for the magnitude in the miss- 
ing bandpass. Positional and magnitude dependent errors 
are seen in the raw photometric data, created by varying 
transmission profile and diffraction effects through the thick 
(5.5 mm) Ha filter, but these are corrected for in the data 
available through the SHS website by comparison with the 
SR data. Photometric consistency is achieved by using the 
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overlap regions between fields to match zero-points across 
the survey. These corrected magnitudes provide a means of 
selecting point-source emitters. The variations in measured 
lAM stellar parameters as a function of field position aris- 
ing from the variable psf from field rotation, vignetting etc, 
especially at large radii from the field centres, requires that 
such selection is performed over limited 1-degree areas. In 
this way stars with an emission line at Ha will show an 
enhanced Ha magnitude compared with the SR magnitude. 
At the bright end of the magnitude distribution, severe pho- 
tographic and SuperCOSMOS saturation effects come in to 
play, limiting stellar photometry to R of about 11-12 in both 
the Ha and SR pass-bands. 

8 SPURIOUS IMAGES IN THE SHS 

Spurious images appear from time to time in all photo- 
graphic images scanned by SuperCOSMOS. They have a va- 
riety of forms and causes and axe present in images extracted 
from the SSS as well as SHS. They can sometimes be picked 
up by examination of the pixel images directly, though they 
are often missed, and can also appear as spurious detections 
in the lAM data. They have a variety of sizes and shapes and 
may be in or out of focus depending on whether the contam- 
inating source is on the emulsion surface or on the platten 
used by SuperCOSMOS to sandwich the film flat for scan- 
ning. Here, we difi^crcntiatc between spurious images in the 
emulsion itself caused by processing defects, emulsion flaws 
and static marks, and those caused by foreign objects on 
the surface of the emulsion or on the back of the fllm. Holes 
and scratches in the emulsion surface can also give rise to 
spurious images. Satellite trails and transient phenomenon 
also give rise to real developed images which may have no 
counterpart in other survey bands of the same region. We 
do not consider these here. 

8.1 Basic causes 

The SuperCOSMOS facility is situated in a class-100 clean 
room and each film is pressure air-cleaned prior to scan- 



ning. However, despite best efforts, particles that may al- 
ready have been present on the emulsion before shipment 
to SuperCOSMOS manifest themselves as spurious images. 
The biggest cause is flne particulate dust (20-100/im). Un- 
fortunately the Estar film base of the Tech- Pan emulsion is 
prone to static charge build-up. 



8.2 Recognising spurious images 

The SuperCOSMOS scanning system is highly specular so 
detritus present on the emulsion surface that is often in- 
visible when viewed under diffuse illumination conditions 
(such as on a light table) is revealed in sharp relief in the 
SuperCOSMOS data. The number of artefacts seen in the 
SHS data is somewhat worse than on other glass plate based 
surveys of the SSS. Fortunately, having matched exposures 
in two bands makes identification of such artefacts more 
straightforward. For example, since the Ha and SR expo- 
sures are registered on the same pixel grid, quotient imaging 
can reveal the locations of spurious images. A 5 x 5 arcminute 
region extracted from h273, a field with a particularly high 
number of spurious images, is shown in Figure 9. 

We can take advantage of the fact that the pixel im- 
age properties of spurious images are usually quite distinct 
from real astronomical images, often having a sharpness be- 
low that possible from the combination of telescope optics 
and seeing disk. Their shapes are often highly irregular and 
non-symmetrical such that they would not fit any normal 
psf. This makes them amenable to Fourier filtering. Objects 
that have no counterpart in the other band are potential 
spurious image candidates though variable objects, novae 
and the effects of de-blending complicate the issue signifi- 
cantly. Various lAM parameters such as the profile statistic, 
ellipticity etc. may also aid in identification. Furthemore, 
spurious lAM objects arising from de-blending overlaying 
contaminating fibres or hairs usually have very high elliptic- 
ities which may help in isolating likely candidates. Storkey 
et al. (2004) discuss techniques for recognising and eliminat- 
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Figure 9. 5 X 5 arcminute extracts of SuperCOSMOS Hq data from SHS survey field h273 highlighting a contaminating spurious image 
together with a matching image with the lAM data overlaid and with all the spurious images in the frame highlighted. 



ing spurious objects in the on-line SuperCOSMOS surveys. 
As yet this procedure has not been apphed to the SHS data. 

9 FLAT-FIELDING OF THE SURVEY DATA 

For any interference filter of the size used here, low-level non- 
uniformities exist which lead to residual non-physical back- 
ground variations in the exposed images. In order to estab- 
lish the magnitude of such effects, three flat-field exposures 
were taken with the filter subject to uniform illumination. 
The flat-field images permitted evaluation of the combined 
effects of filter transmission in the fast, //2.48 Schmidt beam 
and telescope vignetting (see UKST Unit handbook, Trit- 
ton, 1983). The flat-field images were exposed to place them 
on the linear portion of the film's characteristic curve and 
were averaged to give the filter/telescope transmission pro- 
file shown in Figure lOa-b. 

Figure 10a shows transmission contours at 85, 90, 92, 
95, 97 and 98 per cent of maximum transmission in the cen- 
tral region. The response is seen to be asymmetric, with 
the 97 per cent contour extending beyond the edge of the 
5.16° field on the right, which corresponds to the west of the 
survey fields. In the east the transmission decreases more 
rapidly, reaching 85 per cent at the eastern edge. Towards 
the filter corners the transmission drops further. However, 
the 4° overlapping centres (Section 4) and the asymmetric 
nature of the response allows the selection of Ha data re- 
quiring fiat-field correction of less than 15 per cent for any 



given area of sky, provided that adjacent fields are avail- 
able. Most data will require much smaller flat-field correc- 
tions. The effect of fiat-field corrections as large as 15 per 
cent on pixel data is considered in the survey calibration 
section. Figure 10b is a histogram equalization of the actual 
flat-field pixel map which reveals the extremely low level 
artefacts present at the 0.1 per cent level which are invisi- 
ble in a linear rendition. The fiat-field correction has been 
stored as a transmission array with maximum value unity, 
so it is applied by dividing survey image values in inten- 
sity space by the relevant correction array elements. This 
is available as the default option on the SHS website. The 
correction breaks down towards the corners of the scanned 
image and in regions outside of the clear aperture because 
the density of the exposures at the edges of the circular 
aperture is too low to lie on the straight line portion of the 
characteristic curve leading to an over-correction and also 
because the Ha filter transmission is becoming increasingly 
skewed in these extreme regions. Raw transmission or pho- 
tographic density values and generically calibrated intensity 
values without flat-field correction can also be requested on 
the download form. The lAM data is obtained from the raw 
SuperCOSMOS scans without fiat-field correction. 

9.1 Specific filter features 

Despite the superb quality of the filter, low-level, large-scale 
variations in transmitted flux can be seen in the Ha sur- 
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Figure 10. (a - top) Contour plot of the narrowband Ha fil- 
ter transmission over the SuperCOSMOS scanned area of 5.156° 
(1728 X 1728pixel). The contours are plotted at 85, 90, 92, 95, 97 
and 98 jjci cent. The bottom image (b) is a histogram equalisa- 
tion of the same fiat-field data which reveals the filter artefacts 
present at the 0.1 per cent level 

vey images under certain exposure conditions. In particular 
there are tvi^o parallel bands of slightly enhanced transmis- 
sion (leading to elevated photographic density) going El-W 
in the north and south of the filter. These bands are only 
1-3 per cent higher in intensity than the surrounding re- 
gions. A series of low-level artefacts which are not obvious 
in the contour plot because they are at a level of < 1 per 
cent can just be discerned in the filter transmission image 
in Figure lOa-b. They can also just be seen in the contour 
plot as the spike in the 98 per cent contour just right of 
centre. The observed shape of these artefacts mimics the se- 
ries of shallow concentric grooves scored into the surface of 
the mandrel to enable the Tech-Pan film to be sucked under 
light vacuum to the curved focal surface of the plate-holder 




Figure 11. Two 16x blocked down, 5.156° SHS Ha images of 
field h410. The top image (a) is the raw intensity data without fiat 
field correction. The bright areas trace not Galactic emission but 
the filter transmission profile seen in Figure 10. The bottom image 
(b) has had the flat-field correction applied which successfully 
removes the filter artefacts from within the circular aperture. 

to ensure proper focus. They are thought to arise from the 
backscattcring off these grooves of light that has passed right 
through the film. Again, it is gratifying that these artefacts, 
present at the < 0.1 per cent level, are effectively removed 
by application of the fiat-field. 

9.2 Application of the flat-field and correction 
validity 

Field h410, which sits away from the Galactic Plane on the 
extreme edge of the survey at (Z,6) 330.2°, -|-10.28° has 
been chosen to test the vailidty of the fiat-field correction as 
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it contains a very low-level isotropic background of Galactic 
line emission. In the survey image this will bo moderated 
by the filter response. Figure lla-b shows two 16x blocked 
down images of survey field h410. The top image Figure 11a 
is the raw Ha data, before the application of the flat-field 
correction. The structure evident on the field as lighter areas 
is not Galactic emission but matches the filter transmission 
profile. The two horizontal bars are present and the image is 
less exposed towards the corners where the recorded inten- 
sity is lower and the star density also falls. The bright star 
just to the right of centre is e Lupi. 

Application of the flat-field correction results in Fig- 
ure lib which should have a flat background wash of emis- 
sion across it. The structure from the filter is no longer evi- 
dent and the bright areas are now in the corners, where the 
larger flat-field correction over-corrects the SuperCOSMOS 
intensity counts. This will not adversely affect the majority 
of the pixel data available on the SHS website. Data can 
always be taken from the best area of fllter response and 
no flat-fleld correction larger than 15 per cent is necessary 
for any of the survey data that overlap a neighbouring fleld. 
Data from the edge of the survey, where no adjacent fleld 
exists, have been made available and may require a correc- 
tion greater than 15%. Areas affected in this way are flagged 
in the third extension table which accompanies the down- 
loaded FITS image. 

10 GEOCORONAL Ha EMISSION 

Geocoronal Ha emission is caused by fluorescence after solar 
Lyman /? excitation of atomic hydrogen in the exosphere. 
Because imaging surveys lack velocity resolution for the 
emission they record, the geocoronal contribution will be 
present in all of them but indistinguishable from bona fide 
Galactic Ha emission. Fortunately, one modern Ha survey, 
WHAM (e.g. Haffner et al. 2003), offers very good veloc- 
ity resolution (~12 km/s) and is able to separate the atmo- 
spheric emission from the Galactic emission and measure the 
intensity. Nossal et al. (2001) report on Ha observations car- 



ried out by the WHAM instrument in 1997, the same time as 
the SHS imaging was starting at the UKST. They find that 
the geocoronal emission intensity depends on how rrmch the 
earth shades the line of sight from sunlight. Their result- 
ing plot of geocoronal Ha emission as a function of earth 
shadow height shows that at heights greater than 6000 km 
only a very low level ~2 R wash of geocoronal Ha emission 
is present. Based on the observational information available 
in the headers of the SHS images, it is possible to calculate 
the shadow heights for any field. For a random sample of 6 
SHS fields the shadow height for the whole three hour ob- 
servation and across the five-degree field of view was found 
to be greater than 6000 km, so low-level geocoronal emission 
is not problematic, as most Galactic Plane fields covered by 
the SHS will contain significantly stronger emission. 

11 APPLICATION OF AN ABSOLUTE 

CALIBRATION TO THE Ha SURVEY DATA 

The AAO/UKST Ha survey data needs an absolute inten- 
sity calibration if the full scientific value of its sensitivity to 
faint, diffuse emission is to be realized. The intensity calibra- 
tion must provide a reliable means of transforming the pixel 
intensity values from SuperCOSMOS scans of the Ha im- 
ages into meaningful intensity units such as Rayleighs which 
is consistent from field to field. We show that continuum 
emission can be successfully removed from the Ha images 
by scaling and subtracting the SR continuum image. Unlike 
CCD data, which enjoys a linear response over a wide range 
of emission strength, photographic data can be very difHcult 
to calibrate because the response of the emulsion and Su- 
perCOSMOS scanner is linear only over a relatively small 
dynamic range. Variations in sensitivity and background oc- 
cur from exposure to exposure, especially when the Ha and 
SR pairs were taken on different nights, phases of the moon 
etc. Despite this, we show that the survey data have been 
well exposed to capture Galactic emission on the linear part 
of the characteristic curve and can be calibrated by means 
of comparison with the complementary, accurately intensity 
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calibrated SHASSA survey (Gaustad et al. 2001). This pro- 
cess does not form part of the current SHS release but can 
be undertaken by the user as required. 

11.1 Image comparison with SHASSA 

The Southern Ha Sky Survey Atlas by Gaustad et al. (2001) 
provides wide-field narrow-band CCD Ha images of the 
southern sky below 5 = -|-15° taken with a robotic imag- 
ing camera sited at Cerro Tololo Inter-American Observa- 
tory (CTIO). The camera used a small, fast, f/1.6 Canon 
lens which gave a very large (13°) field of view and a spatial 
resolution of 48 arcseconds. Each SHASSA field was imaged 
through a narrow-band interference filter of width 32 A cen- 
tred at 6563 A as well as a continuum filter with two bands 
of 61 A at 6440 A and 6770 A, on either side the Ha line. The 
SHASSA website'* makes available the raw Ha and contin- 
uum images as well as 48 arcsecond and 4 arcminute reso- 
lution continuum subtracted, intensity calibrated data. The 
SHS data is superior in terms of resolution while the general 
sensitivity of both surveys appears qualitatively similar for 
large-scale emission features. For example, in Figure 12 we 
present SHS and SHASSA images of the HII region RCW 19. 

11.2 The SHASSA Intensity calibration 

The SHASSA intensity calibration was derived from the 
planetary nebula spectrophotometric standards of Dopita & 
Hua (1997) after the continuum images had been scaled and 
subtracted from the Ha frames. Aperture photometry for 
eighteen of the bright PNe standards was measured from 
the SHASSA images and used to calculate the calibration 
factor for the whole survey. A difficulty in applying PNe 
line fluxes to Ha narrow-band imaging is the proximity of 
the two [Nil] AA6548, 6584 lines which are included in the 
flanks of the SHASSA Ha filter bandpass. These vary in 
strength relative to Ha between PNe and could significantly 
affect the result. Calculating the transmission properties of 

^ http: / / amundsen.swcirthmore.edu/SHASSA/ 



the interference filter to these linos is complicated by the 
blue-shifting of the bandpass with incident angle, an effect 
which rrmst be treated carefully in the SHS data as the filter 
sits in the fast f/2.48 beam of the UKST. This problem is 
not as severe for the SHASSA data because in this case the 
filter sits in front of the camera lens, leaving only the effects 
of the very large field of view. These effects are considered in 
Section 4 of Gaustad et al. (2001) and carefully accounted 
for in their calibration. 

To allow a more detailed comparison, aperture photom- 
etry for 87 PNe with a range of surface brightness and in- 
tegrated flux and with an independent measure of Ha flux, 
was carried out on SHASSA images. Published spectroscopic 
data were used to deconvolve the contribution from the [N 
II] lines passed by the SHASSA fllter. The results agree 
with pubUshed data to A F(Ha) = -O.Oldex, a = 0.05 for 
SHASSA minus literature fluxes (Prew 2005, in preparation; 
cf. Pierce et al. 2004). Since the PNe literature fluxes have 
associated errors, the SHASSA calibration is better than 
±10 per cent across the whole survey, in agreement with the 
nominal error supplied by Gaustad et al. (2001). 

An additional uncertainty is introduced to the zero- 
point of the SHASSA intensity calibration from the con- 
tribution of geocoronal emission. Gaustad et al. (2001) es- 
timate this by comparison with overlapping WHAM data 
points and interpolating where there are none. Our check of 
the intensity calibration against independent flux measures 
of planetary nebulae indicates that the geocoronal contri- 
bution to the SHASSA Ha images has been successfully re- 
moved. Finkbeiner (2003) also showed there is no signiflcant 
offset between WHAM and SHASSA data. So we conclude 
that the SHASSA data has been well calibrated to a zero- 
point consistent with independent measurements and there- 
fore have confidence in its use as a baseline calibration for 
the SHS data. 




Figure 12. An 105 X 75 arcminute region centred on well Icnown HII region ROW 19 (GUM 10) at a,S = 08''16'"17^, -35°57'58" (J2000) 
extracted from the SHS survey (top left) and the equivalent SHASSA region (top right). A full resolution 8.3 x 4.6 ajrcminute region 
from the SHS survey data is shown below centred on the southern HII region component at 08'' 19"" , -36° 14' and indicated in the lower 
resolution SHS image with a rectangle. 



11.3 Continuum Subtraction of the SHS 

Diffuse emission recorded through the narrow-band Ha filter 
on the Tech-Pan films will be a combination of Galactic Ha 
line emission, continuum emission, night-sky auroral lines 
and geocoronal emission. Ideally all of these components 
would need to be disentangled to extract just the Galactic 
Ha emission. In practice the geocoronal and auroral emission 
is considered as a low (~ 2 R) level but temporally varying 
uniform wash which simply elevates the general background 
on each exposure to a slightly varying degree. 

The matching SR images provide a measure of the con- 
tinuum component and, properly scaled, can be used to pro- 
duce continuum subtracted Ha images. Although the Ha 
and SR exposures are generally exposed to attain the same 
depth for continuum point sources, the nature of photog- 
raphy and the vagaries of the observing conditions (e.g. if 



the exposure pairs wore not contemporaneous and taken in 
different moon phases or if the seeing changed) mean that 
the depth and image quality between the Ha and SR expo- 
sures can and docs vary. Hence it is necessary to determine 
a continuum subtraction scaling factor between Ha and SR 
on a field by field basis. This factor must be precisely de- 
termined for the continuum subtraction to be effective. For 
high-dynamic range CCD exposures, the standard method 
for determining the appropriate scaling factor to subtract 
continuum from narrow-band is to compare aperture pho- 
tometry for stars on both images whose exposures are nor- 
mally interleaved on short timescales. 

Unfortunately, this does not work well with the Ha SR 
film exposure pairs (Pierce 2005), often leading to under or 
over-subtraction of the continuum. This arises due to vary- 
ing backgrounds on the film exposures caused by; low-level 
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emulsion sensitivity variations between films (especially if 
they come from different hypersensitised batches), inherent 
chemical fog variations in the emulsion, processing varia- 
tions and true sky background variations arising for the rea- 
sons given above. These varying backgrounds result in the 
same magnitude stars saturating at different levels on dif- 
ferent exposures as their Gaussian point spread functions 
are superimposed on top of any diffuse emission and ele- 
vated background which can severely truncate their peaks. 
The limited dynamic range of SuperCOSMOS also acts as a 
further low ceiling above the background, which leaves little 
room for these bright stellar Gaussians making it hard to ef- 
fectively utilise stellar photometry to determine the correct 
scaling factors. 

Fortunately, we are able to use the existing SHASSA 
data to provide a well-determined scale with an indepen- 
dently confirmed zero-point to compare with and calibrate 
the SHS images. Even exposure pairs taken years apart can 
be successfully continuum subtracted. A detailed investiga- 
tion of the SHS calibration process has been undertaken by 
Pierce (2005) but the essential aspects of this process and 
its application are given here. For example, Figure 13 shows 
three images of a 30 arcminute region taken from field h350, 
which shows strong, varying Galactic Ha emission. The top 
image is the Ha image downloaded from the SHS website, 
the middle image is its SR counterpart while the bottom im- 
age has had the SR 'continuum' image scaled and subtracted 
via a comparison with SHASSA data. In general this sub- 
traction is very good, removing most of the stellar images 
and the diffuse contiimum. Only the stars which sit in the 
strongest emission have been over-subtracted and appear as 
white spots in the image. 

For a given area of sky, pixel data from each survey 
can be downloaded and, after matching for spatial resolu- 
tion, the SuperCOSMOS intensity counts can be compared 
directly with the Rayleigh values in the SHASSA data. A 
plot of continuum-subtracted SHASSA pixels against equiv- 
alent SHS pixels should return a linear relation with a 
common zero point if the reduction and intensity calibra- 




Figure 13. 30 arcminute region at 16''47'", -49°00' (J2000) from 
near the centre of SHS field h350 with bright areas of emission 
shown darkest. The top image is the raw Ha image, the middle the 
matching SR while the bottom image is a continuum-subtracted 
image of the same area. The subtraction, based on a comparison 
with SHASSA data, has worked well, leaving only minor residuals 
around stars on the bright difi^use emission. 



22 Quentin Parker et al. 

tion have been properly carried out. Comparing incorrectly 
continuum-subtracted SHS data with SHASSA data results 
in an offset between the two surveys. A range of values for 
the scaling factor can be applied to the SHS data until the 
zero-point of the continuum-subtracted UKST survey im- 
ages best matches the zero-point of the equivalent SHASSA 
data, indicating the appropriate value to use. A calibration 
based on the SHASSA data will provide an advantage for 
the SHS over the CTIO survey as it can be applied to the 
full resolution pixel data. This offers the chance to determine 
intensities for emission structures not resolved by SHASSA 
such as the new sample of extended PNe discovered from the 
SHS data (e.g. Parker et al. 2003, 2005, Pierce et al. 2004, 
Prew & Parker, in preparation and see Section 14.1). 

Each scanned SHS survey field, at the full 0.67 arcsec- 
ond resolution, contains over 2 Gb of data so it was not 
practical to download and compare all the pixel data for 
each field. Instead, most of the emission variation on a given 
survey field can be sampled using carefully selected 30 ar- 
cminute regions. For fifteen fields, several 30 arcminute areas 
wore downloaded to sample the complete dynamic range of 
emission present. Once the scaling reliability over a range 
of fiux levels and central aperture locations was established, 
the best single, 30 arcminute region was chosen from the 
whole frame for the rest of the 233 survey fields to provide 
the base calibration for each field. The Ha filter fiat-field cor- 
rection was applied to remove the low-level non-uniformities 
in transmission across the narrow-band filter. 

Data requiring a flat-field correction of up to 15 per cent 
have been shown to be suitable for inclusion in the survey, 
though in most cases data returned from the SHS website 
require less correction than this. For the regions considered 
here, the pixel data only required flat-fleld correction < 3 
per cent in most cases and, in general, no more than 8 per 
cent. For two flelds, h350 and hll09 regions well away from 
the best area of the filter were also selected to confirm that 
the pixel data behaves as expected when larger flat-field cor- 
rections are required. The results are discussed below. 



11.4 The SHS Calibration process 

A detailed description of the SHS survey calibration process 
based on zero-pointing each SHS fleld to SHASSA is given 
in the thesis of Pierce (2005). The essentials of this scheme 
are described here. Each SHS field is completely covered by 
just one 13° SHASSA fleld and data from the best area of 
SHASSA filter response was chosen for comparison with the 
SHS images. For direct comparison, the Ha and SR SHS 
data at 0.67 arcsecond was re-binned to match the 48 arc- 
second SHASSA pixels using the IDL routine Hrebin^. This 
returns the mean value of the 72 x 72 full resolution SHS pix- 
els that constitute a single SHASSA 48 arcsecond pixel, so 
the calibration factor determined from the comparison plot 
applies to the full resolution SHS data. At this coarse reso- 
lution the SR data was scaled and subtracted from the Ha 
imago. Because the correct scaling factor was not yet known, 
a range of scaling factors from 0.4 to 2.0 was applied so that 
the best value could be selected by matching the SHASSA 
zero-point. The equivalent area of SHASSA data was se- 
lected, aligned and trimmed to match using the IDL, Has- 
TROM routine. These images were then compared directly, 
pixel by pixel, with the re-binned, continuum-subtracted 
SHS data to give a plot of SHASSA values in Rayleighs ver- 
sus SuperCOSMOS counts per re-binned pixel. The linear 
portion of the resulting comparison was then fitted to de- 
termine the number of SuperCOSMOS counts per 0.67 arc- 
second pixel per Rayleigh. Averaging SuperCOSMOS data 
in this way only works properly if all the SuperCOSMOS 
pixels are on the linear portion of the characteristic curve at 
the faint end and unsaturated at the bright end. Once the 
SuperCOSMOS dynamic range is exceeded pixel saturation 
arises and the SHASSA to SuperCOSMOS linear relation 
breaks down as the SuperCOSMOS flux becomes increas- 
ingly underestimated (e.g. Phillipps & Parker, 1993) as seen 
in Figure 17. 

The adopted process was followed for the same 30 ar- 
^ Interactive Data Language: http://www.rsinc.com 
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cminute region from SHS field fi350 as sliown in Figure 13. 
The upper image in Figure 14 shows the SHS image with the 
continuum-subtracted at 48 arcsecond resolution. The lower 
image is the trimmed and aligned SHASSA data. Bright 
stars on the continuum subtracted SHASSA images leave 
significant residuals while in the SHS data any stellar resid- 
uals are barely visible and then only for the very bright- 
est stars. The pixel-to-pixel comparison plot for this area 
is shown in Figure 15. Each point on this plot is the pixel 
value from the SuperCOSMOS scan against the Rayleigh 
value from the SHASSA data. Immediately, it can be seen 
that the SHS data follows a tight linear relation with the 
SHASSA values over a range of several hundred Rayleighs. 
The outliers from this distribution at around 4000 Super- 
COSMOS counts correspond to the SHASSA stellar residu- 
als noted from Figure 14. 

A further nine 30 arcminute areas were examined from 
field h350 to cover the whole dynamic range of diffuse emis- 
sion evident in the field. Figure 16 shows the 16 x blocked 
down Ha image for this field with contours of filter response 
overlaid and boxes indicating the areas used. The large area 
of strong emission in the west of this field is caused by the 
UV flux from OB association Ara lA A (Mel'Nik & Efremov 
1995). In all but two cases the areas were selected from the 
best area of the filter response, the exception being the two 
areas extracted from the SE corner which were included to 
test the effect of using pixel data from the edge of the SHS 
survey fields where the flat-field correction is larger. 

Figure 17 shows the resulting SHASSA versus SHS com- 
parison plot with differently shaded points belonging to dif- 
ferent 30 arcminute areas. Almost the full range of the SHS 
data is shown, with the pixel values showing a good linear 
relation to the SHASSA data from the faintest emission on 
the field at ~ 20 Rayleighs right up to ~ 500 Rayleighs. The 
curve in the trend beyond 500 Rayleighs is due to satura- 
tion effects with the SHS data. These points were discarded 
when making the calibration fit. 

The reciprocal of the slope of the linear part of the 
relation provides a calibration factor of 




Figure 14. Two Ha images of a 30' region of SHS field h350 
centred on 16'"47™ - 49° 00" {J2000) as in Figure 13. The top 
image is the SHS image blocked down to match the 48 arcsecond 
resolution of the SHASSA data and continuum subtracted. The 
bottom image is the equivalent SHASSA image, trimmed and 
afigned ready for comparison. Areas of uncancelled bright stars 
are highlighted 

15.1 counts/pixel/Rayleigh to convert the SuperCOSMOS 
intensity counts to Rayleighs for this field. An estimate of 
the error in this calibration is possible, based on the vertical 
scatter of points about this trend ,as a given SuperCOSMOS 
value matches a Rayleigh value in this vertical distribution. 
In this field the Icr scatter is 21 Rayleighs. 

The data taken from the two areas in the SE corner of 
the SHS survey field required fiat-fleld correction of up to 
20 per cent. They behave very well when compared with the 
data requiring less correction, neatly overlaying the main 
trend, with no change in slope or increase in scatter. This 
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Figure 15. The comparison plot of SHS and SHASSA data from 
the area shown in Figure 14. Each point on the plot is a the re- 
binned pixel value from the SHS data against the Rayleigh value 
form the SHASSA data. Note the linear relation and the large 
dynamic range covered by this comparison. The outliers due to 
uncancelled bright stars in the SHASSA data are highlighted. 
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Figure 17. The comparison plot for all ten 30 arcminute areas 
highlighted in Figure 16. All align on the one trend, implying no 
variation in sensitivity across the SHS survey field. The bright 
photographic saturation is seen from ~ 500 Rayleighs. The grey 
points at the faint end of the plot are taken from the area of 
poor filter response in the SE corner of h350 and align well with 
the rest of the data, indicating that, in areas of bright emission, 
flat-field corrections as large as 20 per cent still return reasonable 
intensity values. 




RA (J2I100) 

Figure 16. Blocked down image of survey field h350 with overlaid 
contours at 85, 90, 92, 95 and 97 per cent. The 30 arcminute areas 
used to make the comparison plot shown in Figure 17 are framed 
by the squares which are drawn in black or white to best con- 
trast with their background. The bright HH region to the right is 
NGC 6188, excited by the young open cluster NGC 6193, nucleus 
of Ara OB 1. 

justifies use of the data out to the 15 per cent flat-field con- 
tour where the emulsion records strong emission. 



11.5 SHS sensitivity limit 

While the results from field h350 have shown that the survey 
has been well tuned to the detection of diffuse emission, as 
well as giving an approximate limit to the point at which the 
photographic emulsion saturates, the faint limit has not yet 



been constrained because there is no really faint emission 
on this field. At the faintest point of 20 Rayleighs there is 
no sign of either survey struggling for sensitivity so a field 
exhibiting fainter emission is required to probe the SHS faint 
sensitivity limit. According to the SHASSA data an area of 
the Southern Galactic plane in Monoceros harbours diffuse 
emission that reaches a level as faint as ~ 2 Rayleighs, which 
is ideal to test the faint limit of the SHS data. Figure 18 
shows the 16 times blocked down Ha image of SHS field 
hll09, which covers this area of sky, with contours of filter 
response and the areas selected shown as with Figure 16. On 
this field, nine 30 arcminute regions have been examined as 
two groups: one, labelled 1-5 in Figure 18, from the area 
of best filter response and the other, labelled A-D, probing 
the combined effects of extremely low levels of emission and 
decreased filter transmission. 

The comparison plot from the first group is shown in 
Figure 19. Once again, there is a clearly defined relation 
between the two, although on this faint field the upper sat- 
uration limit is not reached. It is immediately obvious that 
the SHS data can match the SHASSA data right down to 
the faintest level of emission, although the linear response 
of the SHS data is difficult to determine at this low level. 
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Figure 18. Blocked down SHS Ha image of field hll09 with con- 
tours of filter response at 85, 90, 92, 95 and 97 per cent overlaid. 
Two groups of 30' areas are marked by squares and labelled to 
indicate group membership. 



The SHS data are therefore detecting emission structure as 
faint as 2 Rayleighs on this field, although this sensitivity is 
tempered by the scatter evident in the plot and in the exam- 
ination of the areas from the area of poor filter transmission 
discussed below. The linear fit returns a value of 7.7±0.1 Su- 
perCOSMOS counts/pixel/Rayleigh for this particular field 
and provides a reasonable calibration as the Icr scatter to 
the fit is just 6.2 Rayleighs. Note the factor of two difference 
in the slope of the calibration curve for this low emission 
level field hll09 (no emission measure greater than about 
80 Rayleighs seen in the area considered) compared to that 
obtained for high emission field h350 in Figure 17, which 
returned 15.1 counts/pixel/Rayleigh, which is closer to that 
generally obtained for most fields. This serves to emphasise 
the need for individual field calibration due to the variation 
in SuperCOSMOS pixel intensities on a given field, arising 
primarily from variable fog-level, sky background and result- 
ing SuperCOSMOS and emulsion saturation. 

11.6 Fifteen Fields Studied in Depth 

A total of fifteen fields from a wide variety of Galactic en- 
vironments were studied using several 30 arcminute regions 
in each in order to build up a global picture of the survey 
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Figure 19. Calibration plot of the numbered 30 arcminute ar- 
eas from the best area of filter response on field hll09. Points 
from each area are plotted as a different colour. Note the SHS 
data match the SHASSA data down to the faintest intensity at 
~ 2 Rayleighs. The linear fit is also plotted. 



behaviour. The comparison plots for 4 of these are shown in 
Figure 20 with the linear fit overlaid. In each case a clear, 
essentially linear relation can be seen between the SHS and 
SHASSA data. Generally, different 30 arcminute regions fol- 
low a single trend which indicates little variation in emulsion 
sensitivity across the large SHS images. Three of the fields 
examined, hi 75, h350 and h555 show evidence for satura- 
tion at the bright end and fix the bright limit at ~ 500 to 
~ 600 Rayleighs, while none of the fields appear to reach the 
background fog level in the areas that were compared. 

11.7 Calibration of the entire SHS 

From the analysis of the fifteen fields considered above, it 
is clear that each field requires individual calibration and 
that one, well selected, 30 arcminute area of SHS pixel data 
can be compared with SHASSA to give a working estimate 
of the calibration factor and a satisfactory result in most 
cases. On this basis, the 30 arcminute areas of pixel data 
which covered the greatest dynamic range of emission were 
downloaded from each of the 233 SHS fields, re-binned and 
compared with the equivalent area of SHASSA data. The 
results for each field can be found on the SHS web site. For 
each field the position of the area used, the computed scaling 
factor for continuum subtraction, the linear fit and, where 
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Figure 20. Comparison plots for four of the fifteen SHS fields 



appropriate, the coefficients from a third-degree polynomial 
fit are given. The la vertical scatter about the linear fit is 
also quoted to offer an estimate of the error in the calibra- 
tion. Of the 233 survey fields, 76 are relatively featureless 
and exhibit little emission. These are difficult to fit and in 
9 cases the fit failed completely. Forty three of these show 
evidence of the low-level photographic fog. Of the remain- 
ing 157 fields, 122 are well constrained by a linear fit. For 
the other 35 fields the fit can be improved with a low-order 
polynomial relation. Where this is the case the coefficients 
are included in the table on the SHS website. 

11.8 Calibration check of SHS field overlap 
regions 

There is generous overlap between survey fields because 
of the circular aperture of the filter which allows field-to- 
field consistency check of the calibration. For six overlap 
regions between eight fields, a 30 arcminute area was care- 
fully selected from the best possible compromise of filter 
response between two fields, never requiring fiat-field correc- 
tion greater than 15 per cent. The calibration factor calcu- 
lated from the field centre in the best region of filter response 
was applied to these pixel data from the edge. For seven of 
the eight fields the calibration factor determined from the 
linear fit to the data was used. Calibrated and aligned data 
from the two overlapping SHS fields were plotted pixel by 
pixel at 0.67 arcsecond, 10 arcsecond and 48 arcsecond reso- 
lution. If the independently determined calibration applied 
to each field is consistent, the resulting plot of Rayleighs 
from one field against Rayleighs from its neighbour should 
yield a linear relation with slope of unity and no offset. 

From the six fields examined in this way the results 
give good agreement. Four examples are shown in Figure 21 
from comparisons at 48 arcecond resolution. Here the or- 
dinate and abscissa values are in Rayleighs with each axis 
labelled with its field. The fits are quoted on the plots and 
in Table 4. Two comparisons, 459-392 and 459-458, agree 
to better than 10 per cent, three more, 391-329, 458-391 
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and 460-459 also return linear results but agree to just 14, 
29 and 43 per cent respectively. The bottom 3 comparisons 
are for overlapping fields so the good match implies that 
one could construct seamless, large pixel mosaics. For the 
h350 and h349 the linear relation is not so well behaved. 
The slighly curved trend between these two fields probably 
results from the bright saturation evident in the pixel data 
of the overlap region used (note the higher Rayleigh limits 
for this comparison). 



12 COMPARISON WITH IPHAS 

The SHS and IPHAS Ha surveys have areas of over- 
lap at low declinations which permit a direct compari- 
son to be made between the two complementary surveys. 
Figure 22 shows a 3.3 x 2 arcminute region centred on 
18''47'"42.6% -|-01°33'04" which includes the newly discov- 
ered planetary nebula PHR1847-I-0132, taken from a slightly 
shallow SHS survey field hl332 (exposure number HA18088, 
survey grade A2, but exposure time cut short to 168mins 
cf. 180 normally) . The data have been carefully matched in 
terms of co-ordinate projection but not otherwise processed. 
It is clear the two surveys achieve similar depth for diffuse 
emission but that the IPHAS survey goes deeper for point- 
sources due to its better resolution. Further details of the 
IPHAS survey are given by Drew et al. (2005). 



13 COLOUR-COLOUR PHOTOMETRIC 
PLOTS 

The SR, Ha 'R-equivalent' and I magnitudes can be com- 
bined to provide Ha — R and R — I colours for the objects 
detected by the lAM software. The narrowband Ha pho- 
tometry should be sensitive to point-source emitters, so the 
SHS stellar photometry is of particular interest. A colour- 
magnitude diagram (CMD) of SR magnitude versus Ha — R 
can be constructed from the survey data to trace the average 
values of Ha — R for normal stars by brightness and help 
identify emitters by highlighting objects which stand apart 
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Comparison of fields HAL0459 and HAL0392 
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Comparison of fields HAL0459 and HAL0458 
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Comparison of fields HAL0460 and HAL0459 
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Figure 21. Comparisons between adjacent fields which have been 
calibrated according to a linear fit with SHASSA data. 
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Table 4. Fits obtained from the comparison of SuperCOSMOS pixel data in overlap regions between adjacent survey fields using linear 
calibration factors 



Fields 


Position 


Slope 


SlopelO 


Slope48 


Calibl/Calib2 


350 and 349 


16 38 


-48 


GO 


0.88" 


0.95" 


0.97" 


1.03 


391 and 329 


08 25 


-46 


30 


1.15 


1.13 


1.14 


0.62 


458 and 391 


08 22 


-42 


15 


0.49 


0.54 


0.57 


1.44 


459 and 392 


08 48 


-42 


15 


0.86 


0.90 


0.93 


1.17 


459 and 458 


08 36 


-40 


30 


0.90 


0.91 


0.89 


0.93 


460 and 459 


08 56 


-41 


15 


0.60 


0.71 


0.71 


1.22 



" This result is from a linear fit whereas a polynomial might be more appropriate. 
Refer to top panel in comparison plots given in Figure 21 
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Figure 22. A 3.3 X 2 arcminutc comparison region be- 
tween SHS (top) and IPHAS (bottom) data centred on 
18^47"42.6^+01°33'04" including the new planetary nebula 
PHR1847+0132. 



from this. Figure 23a shows a CMD constructed from a 1° 
region from the centre of survey field hl35. Most stars can be 
seen to congregate around a stellar locus running vertically 
in the diagram at Ha — ~ —0.2 and with an increasing 
spread towards fainter magnitudes. This is due to the census 
including a more complete population of objects at fainter 
magnitudes and increasing photon errors. Outliers can be 



seen either side of the distribution and, if present, emission 
line objects will sit to the left of the main stellar locus. 

A colour-colour plot of i? — / vs Ha — R provides more 
information about the stellar population along a given line 
of sight. Figure 23b presents such a plot for the same area 
of sky as described above and uses the / band magnitudes 
which have been corrected to the SR photometry. In this 
plot the stellar locus is centred at Ha — ~ —0.2 and 
R — I =0.15. The area covered suffers high reddening of 
up to E{B — V)= 4.7 mag (Schlegel, Finkbeiner, & Davis, 
1998) and this is evident in the stretching of the stellar lo- 
cus towards larger values of R — I and, significantly. Ha — R 
making them appear to be emission line stars. These red- 
dened non-emitters can be identified in the colour-colour 
plot of R — I versus Ha — R and excluded from studies of 
potential point-source emitters (Pierce 2005). An additional 
complication is the potential contamination from late type 
stars. A TiO opacity minimum near 6536A enables the con- 
tinuum to be attained, producing the peak, compared to 
the TiO band heads on either side of the Ha filter. Fur- 
ther to the red, the Tech-Pan emulsion sensitivity cuts off 
at 6990A. Such objects can thus appear as apparent Ha 
emitters when compared to the matching SR photometry 
unless the complementary I band photometry is included, 
as such late type stars will be brighter in this band than Ha 
emitters. However, colour-colour plots created in 1-degree 
sub-regions, avoids the smearing effects on the photometry 
due to small positional shifts in the stellar locus across a Ha 
survey field. These have proved very effective in identifying 
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Colour magnitude diagram for HAL0135 
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Figure 23. The colour magnitude diagram of SR vs Ha — R and 
the colour-colour plot of — / vs Ha — R for stars from an area 
centred on the middle of SHS field hl35. On average the stars 
sit at roughly neutral values of Ha — R = -0.2 and R — 7 = 0.15. 

The reddening in this field of up to E(B-V) = 4.7 mag is evident 
in the colour-colour plot as a smearing of the stellar locus towards 
higher values of R — I and Ha — R. The circled point indicates 
the location of a Planetary Nebulae. 



point-source emission candidates and has been successfully 
employed to provide targets for follow-up multi-object spec- 
troscopy with 6dF and 2dF at the AAO (e.g. Hopewell et 
al. 2005). 



nebula project (MASH; Parker et al. 2003, 2005 in prepa- 
ration) which has uncovered about 1000 new Galactic plan- 
etary nebulae (PNo), nearly doubling the sample accrued 
from all sources over the last 100 years. Related projects 
concern identification of a significant new PNe population 
in the Galactic Bulge (e.g. Peyaud, Parker & Acker 2003) 
and the discovery of an important sample of Wolf-Rayet 
central stars of PNe (e.g. Morgan, Parker & Russeil 2001; 
Parker & Morgan 2003) including the detection of the only 
[WN] PN central star in the galaxy (Morgan, Parker & Go- 
hen 2003). A possible new phase of PNe evolution has also 
been reported around a strongly masing OH-IR stax (Cohen, 
Parker & Chapman 2005). A very large PN in an early stage 
of interaction with the ISM has also been discovered (Pierce 
et al. 2004) as well as two, very large bipolar PNe previously 
mis- identified as HH regions (Prew, Parker & Russeil 2005). 
A new sample of large (>4 arcminute), highly evolved, low 
surface brightness PNe have also been found from exami- 
nation of 16 X blocked down FITS images of the entire 233 
fields of SHS survey. These blocked-down images effectively 
enhance large angular size low surface brightness features 
(e.g. Pierce et al. 2004; Frew & Parker 2005, in preparation). 
All these discoveries are being investigated with follow-up 
spectroscopy to determine the fundamental parameters of 
this significant new sample. 



14 SCIENTIFIC EXPLOITATION 

The SHS on-line atlas was released in stages starting in 2002 
with the complete survey made available in 2003. A variety 
of programmes are underway to exploit the scientific poten- 
tial of this new resource. Several illustrative project exam- 
ples are briefiy mentioned below. 



14.1 Planetary Nebulae 

The largest project arising from the AAO/UKST Ha sur- 
vey has been the Macquarie/AAO/Strasbourg Ha planetary 



14.2 HII regions and regions of star formation 

This type of study is currently ripe for exploitation with lit- 
tle work currently undertaken. Mader et al. (1999) report the 
discovery of a significant new population of Herbig-Haro ob- 
jects from an SHS extension field in Orion and a large wind- 
blown bubble with secondary star-forming regions around 
its periphery by Cohen et al. (2002). Numerous very faint 
HII regions have also been discovered as a by-product of the 
MASH survey (Parker et al. in preparation). Many large 
features are evident including bubbles which trigger star- 
formation and induce velocity departures of the associated 
HII regions. Using data from the SHS and additional kine- 
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matic information from Fabry- Perot observations (Georgelin 
et al. 2000), the filamentary Ha counterparts and triggered 
HIl regions for the HI shell centered at 290.1+0.2 (Rizzo 
and Arnal, 1998) have been revealed. The high resolution 
of the data permits precise description of the morphology 
and extent of such HH regions. This is vital information 
in determining the location of the exciting stars which can 
be inferred via orientation of the observed rims and dust 
"elephant trunk" with respect to the HII region as a whole. 
Furthermore, the Ha counterpart and visible extension can 
be directly compared to radio HH regions. This is essential 
information needed to determine the distance of HH regions 
in the framework of the study of the large scale structure of 
our Galaxy (Russeil et al., 2005). 

14.3 Supernova Remnants 

Several programmes searching for the optical counterparts of 
supernova remnants (SNRs) in the SHS have already been 
undertaken. Walker, Zealey & Parker (2001) report find- 
ing new filamentary shell structures traced by Ha emission 
that arc likely associated with Galactic SNRs. A more recent 
project is underway to uncover SNR candidates across the 
entire SHS from careful scrutiny of both the blocked-down 
FITS images and the original survey films, with several new 
SNRs already confirmed (Stupar, private communication). 
Searches for new optical Ha counterparts around the known 
Galactic SNR overlapping the SHS is also underway. One 
new Galactic SNR discovered serendipitously via the MASH 
programme has already been reported (Parker, Frew & Stu- 
par 2004) . A significant increase in the known population of 
optically detected Galactic SNRs is promised. 

14.4 Point-source emitters 

One area of more recent study is the search for point-source 
emitters and the subsequent follow-up spectroscopy of can- 
didates identified from the SHS Ha and SR photometry. 
Drew et al (2004) report the discovery, via SHS photome- 
try, of only the 4*^ known massive WO star in the Milky 



Way Galaxy identified as part of a general programme of 
candidate point-source follow-up. Additionally, Hopewell et 
al. (2005) present five new WC9 stars discovered from the 
SHS data in a similar fashion. Pierce (2005) and Pierce et 
al. (2005, in preparation) demonstrate the power of the SHS 
to reveal significant new populations of Ha emitters via a 
particular study in the Vela molecular ridge, especially when 
combined with I band and 2MASS photometry. These pre- 
liminary projects have been finding, for the magnitude range 
explored most thoroughly (12 < R < 16.5), that 10-20 per 
cent of candidates are confirmed as emission line objects via 
follow-up spectroscopy. Their Ha equivalent widths usually 
exceed 20A. More recent work by Pierce (2005) indicates 
that the situation can be improved by weeding out M stars 
more thoroughly using the 2MASS data. 

15 CONCLUSIONS 

The AAO/UKST Ha survey as scanned by SuperCOSMOS 
is now complete and on-line as the SHS atlas. It represents 
a powerful tool for the study of the ionized gas content 
of our galaxy on a range of spatial scales from arcsecond 
to tens of degrees. The distribution and structure of the 
ionised gas result from a wide range of astrophysically inter- 
esting phenomena. The astrometric and photometric prop- 
erties have been described and shown to be well behaved 
and adequate for most purposes. Importantly, despite dif- 
ficulties associated with photographic data and the scan- 
ning process, comparison with the independently calibrated 
SHASSA images has shown that the SHS survey faithfully 
records diffuse Galactic emission over a wide range of inten- 
sities from ~ 5 Rayleighs to 500 Rayleighs. Emission down to 
~2Rayleighs has been detected on one field, hi 109. A cali- 
bration scheme for all 233 survey fields has been generated, 
based on comparison of a carefully selected, 30 arcminute re- 
gion from each field with the equivalent area of intensity cal- 
ibrated SHASSA Ha image. If the limitations of the data are 
respected in terms of dynamic range, reliable flux estimates 
are possible. The survey is clearly appropriate for studies 
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of individual Ha emitting objects including point-sources as 
well as being suitable for the study of the ionized interstel- 
lar medium in general. A variety of projects exploiting this 
resource are already underway and many exciting discover- 
ies have already been made. The community is invited to 
consider use of this valuable survey when undertaking any 
study of the Southern Galactic Plane. 
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